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Abstract 

A method is developed for dealing with ultraviolet divergences in calcula- 
tions of cosmological correlations, which does not depend on dimensional 
regularization. An extended version of the WKB approximation is used to 
analyze the divergences in these calculations, and these divergences are con- 
trolled by the introduction of Pauli-Villars regulator fields. This approach 
is illustrated in the theory of a scalar field with arbitrary self-interactions 
in a fixed flat-space Robertson-Walker metric with arbitrary scale factor 
a(t). Explicit formulas are given for the counterterms needed to cancel all 
dependence on the regulator properties, and an explicit prescription is given 
for calculating finite regulator-independent correlation functions. The pos- 
sibility of infrared divergences in this theory is briefly considered. 
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I. INTRODUCTION 



Much effort has been expended in recent years in the calculation of quan- 
tum effects on cosmological correlations produced during inflation. These 
calculations are complicated by the occurrence of ultraviolet divergences, 

which have typically been treated by the method of dimensional regular- 
ization. Unfortunately, this method has several drawbacks. It is difficult 
or impossible to employ dimensional regularization unless the analytic form 
of the integrand as a function of wave number is explicitly known, so cal- 
culations have generally relied on an assumption of slow roll inflation, or 
even strictly exponential inflation. Also, even where an analytic form of the 
integrand is known, dimensional regularization can be tricky. Senatore and 
Zaldarriaga[l] have shown that there are terms in correlation functions that 
were omitted in work by other authors [2], [3]. 

This article will describe a method of dealing with ultraviolet diver- 
gences in cosmological correlations, without dimensional regularization. For 
the purposes of regularization of infinities, we employ a generally covariant 
version of Pauli-Villars regularization [4]. In order to calculate the coun- 
terterms that are needed to cancel infinities when the regulator masses go 
to infinity, we introduce an extended version of the WKB approximation 
(keeping not only terms of leading order in wavelength), which works well 
even when the wave number dependence of the integrand is not explicitly 
known, and can therefore be applied for an arbitrary history of expansion 
during infiation. 

This method is described here in a classic model, the fluctuations of a 
real scalar field in a fixed general Robertson- Walker metric. This is simple 
enough to illustrate the use of the method without the general idea being 
lost in the complications of quantum gravity, and yet sufficiently general 
so that we can see how to deal with an arbitrary expansion history. As 
we shall see, these methods yield a prescription for calculating correlation 
functions that are not only free of ultraviolet divergences, but independent 
of the properties of the regulator fields. 



II. THE MODEL 



We consider the theory of a single real scalar field ip{x) in a fixed metric 
Qfivix), with Lagrangian density 



L = ^J-T)eig 



(1) 
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where V{ip) is a general potential. The modifications in this Lagrangian 
needed to introduce counterterms and regulator fields will be discussed in 
Sections III and IV, respectively. 

This theory will be studied in the case of a general flat-space Robertson- 
Walker metric: 

500 = -1 , ffoi = , gij = a^it) dij , (2) 

with a{t) a fixed function (unrelated to V{ip)), which is arbitrary except 
that we assume that a{t) increases monotonically from a value that vanishes 
for t — > — GO. The field equation is then 

(f + 3Hip - a-^VV + y'iv) = , (3) 

where as usual H = a/a is the expansion rate. We define a fluctuation S(p 
by writing 

ip{^,t) = ^{t) + 6ip{x,t) , (4) 
where (f{t) is a position-independent c-number solution of the field equation: 

^ + mip + V'{Tp) = . (5) 

Our calculations will be done using an interaction picture, in which the 
time-dependence of 5tp is governed by the part of the Hamiltonian quadratic 
in 5ip, so that dip satisfies a linear differential equation 

5(p + m5^-a-'^V'^8^ + V"{Tp)5ip = Q . (6) 

The commutation relations of S^p are 

[5^{-^,t),5<p{y,t)]=ia-\t)5\-K-y) , (7) 

[<5<p(x, t),5^{y, t)] = [(5vj(x, t),6^{y, i)] = . (8) 
The fluctuation can therefore be expressed as 

where Q;(q) is an operator satisfying the familiar commutation relations 

[a(q), «t(q')] = 5\ci - q') , [a(q), «(q')] = , (10) 
and Uq{t) satisfies the differential equation 

Uq + 3HUq + a'^q^Uq + V"ilp)Uq = (11) 
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and the initial condition, that for t — >• — oo, 



(2^)3/2a(t)V2^ 



exp 



tq 



^ dt' 



a{t') 



(12) 



where T is an arbitrary fixed time. (The commutation relations (10) and 
the initial condition (12) ensure that the commutation relations (7) and 
(8) are satisfied for f — t- — oo. The three commutators in these commutation 
relations satisfy coupled first-order differential equations in time, which with 
this initial condition imply that the commutation relations are satisfied for 
all times.) 

According to the "in-in" formalism[5], the vacuum expectation value of 
a product Onit) of Heisenberg picture fields and their derivatives, all at 
time t, is given 



{Onit)) 



VAC 



Texp 



H'j{t')dt')Oi{t) Texp 



H'j{t')dt' 



(13) 

where (• • •)o denotes the expectation value in a bare vacuum state annihi- 
lated by a(q); T and T denote time-ordered and anti-time-ordered prod- 
ucts; Oj{t) is the operator 0{t) expressed in terms of interaction picture 
fluctuations; and H'^ is the interaction Hamiltonian, the sum of terms in 
the Hamiltonian of third and higher order in the fluctuations, expressed in 
terms of the interaction-picture fluctuation 5ip: 



d-^x 



(14) 



We will evaluate Eq. (13) as an expansion in the number of loops. If 
we like, we can introduce a loop-counting parameter g by writing V{(p) = 
g~'^F{gip), with F{z) a (^-independent function of z, so that the number of 
factors of fl' in a diagram with L loops and E external scalar lines is 



# = 2L-2 + E . 



(15) 



Thus an expansion in the number of loops is the same as a series in powers 
ofg'. 



^It will be implicitly understood that the contours of integration over time are distorted 
at very early times to provide exponential convergence factors, as described in ref. [3]. 
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III. ONE-LOOP COUNTERTERMS 



Infinities are encountered when calculating loop contributions to (13) in 
this model. As in flat space, they can be canceled by introducing suitable 
counterterms into the Lagrangian. (When regulator fields are introduced, 
the counterterms instead cancel dependence on the regulator properties.) 
But the Lagrangian cannot know what metric will be adopted, or the clas- 
sical field ^ around which the field is to be expanded, so neither can the 
counterterms. Thus we must return to the generally covariant form (1) of 
the Lagrangian in analyzing the possible counterterms that may be needed 
and employed. 

The general one-loop one-particle-irreducible diagram consists of a loop 
into which are inserted a number of vertices, to each of which is attached any 
number of external lines. An insertion with external lines is given by the 
(A^-|-2)th derivative of T^(v') with respect to (/? at 93 = ^, so the counterterm 
in the Lagrangian can only be a function of !/"((/?), and of gp,y and its 
derivatives. Furthermore, the operators appearing in a counterterm needed 
to cancel infinities can only be of dimensionality (in powers of energy) four 
or less. But V" (ip) has dimensionality two, so the only generally covariant 
counterterm satisfying these conditions is of the forrro 



(16) 



where R is the usual scalar curvature, and A, B, and C are constants that 
depend on the cutoff (that is, on the regulator masses), but not on the 
potential. Dimensional analysis tells us that in the absence of regulator fields 
A is quadratically divergent, while B and C are logarithmically divergent. 

If we now specialize to the Robertson- Walker metric (2), and write the 
scalar field as in (4) , this counterterm becomes (aside from a c- number term) 



£1 loop 



A{v"'iTp)6^ + ^V""m5^' + ...) 
+B {2V"{ip)V"'m6^ + \V"'''m + V"{ip)V""m6^'' + ...) 



-{6H + 12H^)C {v"'{-^)5^ + -y""((^) V 



+ 



(17) 



This argument does not rule out an additional term proportional to 



\/— Detg V" {ip)g^^ d^ipdvf, but one-loop diagrams do not generate ultraviolet divergent 
terms with spacetime derivatives acting on external line wave functions. 
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These terms are of one-loop order, and hence to that order are to be used 
only in the tree approximation, with a new term in the interaction Hamil- 
tonian given by 

AHi = -Jd^x /:^^°°P . (18) 

The terms shown explicitly in Eq. (17) are the only counterterms in Eq. (18) 
that contribute in one-loop order to the one-point and two-point functions. 

IV. REGULATORS 

The counterterm (17) is certainly not the most general counterterm that 
would be consistent with the symmetries of the Robertson-Walker metric. 
For instance, if we didn't know anything about general covariance, we would 
have no reason to expect that H and should occur in the linear com- 
bination R = —6H — 12H^. In order to be sure that the divergences we 
encounter will be of a form that can be canceled by the counterterm (17), 
although we do our calculations for the Robertson-Walker metric (2), we 
shall adopt a regulator scheme derived from a generally covariant theory. 

The usual approach to this problem is to use dimensional regularization, 
which we wish to avoid for reasons given in Section I. There are other meth- 
ods of regularization that have been extensively applied to the evaluation of 
expectation values of operators like the energy-momentum tensor in curved 
spacetimes[6] but not as far as I know to the calculation of cosmological 
correlations. 

One such method is covariant point-splitting [7]. This method is well 

suited to the calculation of expectation values of bilinear operators, where 
the ultraviolet divergence arises from the confluence of the arguments of the 
two operators. Because it is a covariant method, it can be implemented by a 
renormalization of the bilinear operator that respects its transformation and 
convergence properties. It seems difficult to apply covariant point-splitting 
to the calculation of cosmological correlations, where one integrates over the 
separation of the spacetime arguments of the interaction Hamiltonian.' 

There is another widely used method known as adiabatic regulariza- 
tion[8]. In this method, one subtracts from the integrand its asymptotic form 
for large wave numbers, as determined by an extended version of the WKB 
method. Experience has shown that though not covariant, this method 
yields the same results for expectation values of bilinear operators as covari- 
ant point-splitting [9]. But adiabatic regularization affects the contribution 
of small as well as large internal wave numbers, so it seems unlikely that 
it can be applied to the calculation of cosmological correlations, where for 
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some diagrams the contribution of small virtual wave numbers to correla- 
tion functions depends in a complicated way on external wave numbers, so 
that adiabatic regularization cannot be implemented by the introduction of 
generally covariant counterterms in the Lagrangian. 

We will instead here employ a generally covariant version of Pauli-Villars 
regularization [4], which like covariant point splitting and adiabatic regular- 
ization has previously been applied to the calculation of expectation values. 
For the theory studied here, the Lagrangian (1) is modified to read 

n 

(19) 

where Xn are regulator fields, and Z„ and M„ are real non-zero parameters. 
In order to eliminate ultraviolet divergences up to some even order D, we 
must take the Z„ and regulator masses M„ to satisfy 

n n n n 

(20) 

For instance, if there were only logarithmic divergences then D = 0, and we 
would only need one regulator field, with Zi = —1. In one-loop calculations 
the maximum degree of divergence is quadratic, i.e. D = 2, and to satisfy 
the conditions (20) we need at least two regulator fields. In our calculations 
we will not need to make a specific choice of the number of regulator fields, 
but only assume that there are enough to satisfy Eq. (20). 

The coefficients A, B, and C in the one-loop counterterm (17) will be 
given values depending on the Z„ and M„, such that all expectation values 
(13) approach finite limits independent of the Z„ and Mn, as the M„ become 
infinite. As we will see, this condition not only fixes the terms in A, B, and 
C that are proportional to logarithms of regulator masses and the term in 
A that is proportional to squares of regulator masses, but also the terms 
in A, B, and C that depend on regulator masses only through their ratios, 
and hence that remain fixed as the regulator mass scale goes to infinity. 
The only terms in ^, B, and C that will not be fixed by this condition are 
finite terms independent of regulator properties, which of course represent 
the freedom we have to change the parameters in the potential or to add a 
non-minimal coupling of the scalar field to the curvature. 

The regulator fields Xn like the physical field (p are written as classical 



C = ^/-Detg 
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fields plus fluctuations 

Xn{^,t)=Xn{t) + SXn{^,t) . (21) 

The classical fields satisfy the coupled field equations 

iP + 3Hip + V'(if + J2xn)=0 (22) 

n 

% + ZHYn + Z-^V'{tp + Y, Xn) + Mlxn = ■ (23) 

n 

We assume throughout that the regulator masses M„ are all much larger 
than H(t') and V" (Tp{t')] over the whole range from t' — t- — oo to the 



time t' = t at which the correlations are measured. In consequence, the 
classical field equations (22) and (23) have a solution in which all the Xn 
are less than Tp by factors of order H"^ /M"^ and \V" {jp)\/M'^, and so may be 
neglected. adopt this solution for the classical fields. In particular, the 
field Tp then satisfies the original classical field equation (5) . 

In dealing with internal lines, it is convenient to lump together the phys- 
ical field fiuctuation 5(p and the fluctuations 5xn in the regulator fields, by 
introducing an index N (and likewise M, etc.) such that 8xn is the phys- 
ical field fluctuation 5(p ioi N = Q and is a regulator field fiuctuation for 
N = n > 1, both in the interaction picture. The general field fluctuations 
satisfy the coupled field equations 

SxN + SH6xN -a-^VHxN + MjjdxN + Z],^ V" (^) ^ ^xm = , (24) 

M 

where Zq = 1 and Mq = 0. The commutation relations of the Sx are 

[6xN{^,t),6xM{y,t)] = ia-''{t)6^{x - y)Z^^(5iVM , (25) 

[SxN (x, t) , SxM (y , t)] = [SxN (x, t) , SxM {y,t)] = o. (26) 

The general fluctuation can therefore be expressed as 

5xNi^, «) = E / d'l [e*'i-"aM(qX,(i) + e-*'i-ai,(q)n|J^ (t)] , (27) 

where Q;iv(q) satisfy the commutation relations 

[aiv(q), aL(q')] = -^'(q - qO^^'-^iVM , [ajv(q), aM(q')] = , (28) 
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and the u^g{t) are solutions of Eq. (24): 



L 

distinguished by the initial condition, that for t — )■ — oo, 



(29) 



1 



where 



(27r)3/2a3/2(t)^2«iV5(i) 



KNq{t') = 



-I I KNq{t')dt' 



1/2 



(30) 



(31) 



The aiv(q) are all taken to annihilate the vacuum. The two-point functions 
appearing in propagators are then given by 



(5Xiv(xi,ii)5xM(x2,t2))o = E/ e"i<^^-^'^Z^'u%^{h 

K •' 



U 



K* 
Mq 



(32) 

In calculating one-loop graphs, we must integrate over one or more times 
ti associated with vertices, and over a single co-moving wave number q. 
There are two ranges of q = |q| where the integrand is greatly simplified. 

In the first range, q/a{t) (and hence all q/a{ti)) is much greater than 

1 /2 

H{t') and V" {Tp{t')^ for all t' < t, as well as much greater than the 
physical wave numbers associated with external lines, though q/a{t) is not 
necessarily greater than the regulator masses. In this range, we can reli- 
ably evaluate the integrand in an extended version of the WKB approxi- 
mation, described in an Appendix. Any term that would be convergent in 
the absence of cancelations among the physical and regulator fields makes a 
negligible contribution to the integral over this range. 

In the second range, q/a(t) is much less than the regulator masses, 

though it is not necessarily less than H{t') or V" {Tp{t')^ ^ or the physical 
wave numbers associated with external lines. In this range, it is safe to ig- 
nore the regulator fields. (We do not have to worry about the contribution 
of times t' so much earlier than t that q/a{t') is of the order of the regulator 
masses, because this contribution is exponentially suppressed by the rapid 
oscillation of the integrand at these early times.) 

It is crucially important to our method of calculation that, because we as- 

sume that the regulator masses are much larger than H{t') and V ' 
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and the physical wave numbers associated with external lines, these ranges 
of wave number overlap. We can therefore separate the range of integra- 
tion of co-moving wave number by introducing a quantity Q in the overlap 
region, so that Q/a{t) is much less than all regulator masses, and much 

greater than H{t') and V"(Tp{t')j and the physical wave numbers asso- 
ciated with external lines. We can evaluate the integral over q < Q ignoring 
the regulators, and over q> Q hj using the WKB approximation. No errors 

are introduced by this procedure in the final result, because we are taking 
the regulator masses to be arbitrarily large compared with Q/a{t), which 

1 /2 

is taken to be arbitrarily large compared with H{t') or V" (Tp{t')j for 
t' <t oi the physical wave numbers associated with external lines, so terms 
proportional to quantities like Q/Mna{t) or Ha{t)/Q are entirely negligible. 

It should be emphasized that Q is neither an infrared nor an ultraviolet 
cutoff, but simply a more-or-less arbitrary point at which we choose to split 
the range of integration. As long as Q is chosen in the overlap of the two 
regions defined in the previous paragraphs, the sum of the integrals over 
q < Q and q> Q will automatically be independent of Q. 

V. THE TWO-POINT FUNCTION 

To demonstrate the use of the methods described in the previous section, 
and to evaluate the coefficients A, B, and C in the counterterm (16), we 
will now calculate the one-loop corrections to the vacuum expectation value 
of the product Sipniy, t) Sifniz, t) of Heisenberg picture fields. This can be 
written in terms of a Green's function Gp{t), as 

{Sfniy, t) (z, t))vAc = j '^^P (^P " (y - ^)) ^pi^) ■ (33) 

Leaving aside vacuum fluctuations and countertcrms, there are three one- 
loop diagrams, shown in Figure 1. In this section we will consider only 
the one-particle-irreducible diagrams, I and II; these will suffice to allow us 
in Section VI to fix the coefficients A, B, and C in the counterterm (16). 
Diagram III will be dealt with in Section VII. 

Diagram I 

By the usual rules of the "in-in" formalism, after integrating over spatial 
coordinates, the contribution of diagram I to the two-point function is 

G^(i) = -2(27r)6Re £ dtia^{ti)V"'(ifiti)) (it2 a3(t2) ^"(^(ts)) 
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Ill 



Figure 1: Diagrams for the two-point function. 



E 



7^17-1 



KLMNM'N' 



'"p(*)l^'"p(*l)Wp(*2)«Mq(*l)^M'q(*2)njf^'(*l)^Ar'9'(*2) 



K 



,L* 



, (34) 



where q = |q| and (7' = |q — p|- The first term in the square brackets arises 
from diagrams in which the vertices come either both from the time-ordered 
product or both from the anti-time-ordered product in Eq. (13), while the 
second term arises from diagrams in whicli one vertex comes from the time- 
ordered product and the other from the anti-timc-ordcrcd product. 

As described at the end of the previous section, to calculate Gp{t) we 
divide the region of integration over g = |q| into the ranges q < Q and 
q > Q, where Q is chosen so that Q/a{t) is much less than all regulator 
masses but much greater than p/a(t) and H(t') and \V" {lp{t'))\^/'^ for all 
t' < t. 

For q < Q, we can ignore the regulators, and set K, L, M, N, M', iV' 
all equal to zero, with Uq^ just equal to the wave function Uq in the absence 
of regulators. This contribution takes the form 

Gl'<^(t) = -2(27r)6Re j' dh (h) V" j' dt2 a^{t2)V"' {Tp{t2 



e{ti - t2)ul{t)u;{h)u;{t2 



q<Q 



d q Uq{ti)u*q{t2)Uq'{ti)u*q,{t2) 



-^\Up{t)\'^U*{ti)Up{t2) J d^q U*{ti)Uq{t2)u*,{ti)Uq>{t2) 



(35) 



No limit has been put on the second integral over q, because the oscillat- 
ing exponentials in Uq and Uqi make this integral converge[3], so that the 

contribution of wave numbers with q > Q is exponentially small. 

For q > Q, we can use the WKB approximation (30). This contribution 
then takes the form 

G^'>«(t) = -'2Re dii V^'"(^(ti)) y^*' di2^'"(^(t2))iip(t)<(ii)<(t2) 

d^q 



) I^Kg {t2 ) KLq (* 1 ) I^Lg {t2 ) 



X exp ( -i / [KKg{t') + HLq{t')] dt' ) . 



t2 



(36) 
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Note that we have dropped the distinction between q' and g, because p is 
negligible compared with q for q > Q. We have also dropped the contribution 
of the second term in Eq. (34), because this term converges for each K and 
L, and so makes a negligible contribution to the integral over values q > Q. 

The contribution of values of t2 at any fixed time less than ti is also 
negligible, because of the rapid oscillation of the final factor. But there 
is an important contribution from values of t2 that are so close to ti that 
(^1 — t2)Q/a{ti) is not large. This contribution can be evaluated by setting 
t2 = tl everywhere except in the range of integration in the exponential, so 
that 

Gl'>Q{t) = -2Re f_^dhV"'(ipiti)) £ dt2V"'(if{h))ulit)u;\h) 



' —OO 

dh 



X 

KL 



lq>Q 4=KKq{tl)KLq{tl) 
X exp - t2)[KKqitl) + KLqih)]) 

= - f dtiV"'(lp{h)fl^\ul{t)uf{ti 

J — oo V / L 

\^ Z'^Z'^ [ — (37) 

^ ^ ^ Ja>Q'il^Ka{tl)nLa{tl)\KKa{tl) + t^La{tl)\'^ 



X 

KL 



The integral over q converges because ^fc^ — 0- This integral receives 
contributions from terms where XK and xl are both regulator fields Xm 
and Xni or are a regulator field Xn and a physical field xo = V'; or are two 
physical fields. Adding these contributions gives 

G^'>«(t) = TT j'_^dti a\ti)V"' {Tp{ti)f [ul{t)uf{h) 

l^^m^n M2-M2 +^2^^n ^"^^- + ^"1^ 

Note that, because J2n — this is independent of the units used to 
measure Q and the regulator masses, as long as the same units are used in 
all logarithms. 

Diagram II 

By the usual rules of the "in-in" formalism, after integrating over spatial 
coordinates, the contribution of diagram II to the two-point function (33) is 
given by 

G^^(t) = (27r)3|_*^dtia3(ti)y""(^(ti))lm(«2(^^^*2(^^^^ 
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X J2 Z^' d'qu§^{h)u§Tg{h). (39) 
KNN' ^ 

We again divide the range of integration over g = |q| into the ranges q < Q 
and q > Q, where Q is chosen so that Q/a{t) is much less than ah regulator 
masses but much greater than p/a{t) and H{t') and \V" {f{t'))\^/'^ for all 
t' < h. 

For q < Qwe can ignore the regulators, and set K, N, and N' all equal to 
zero, with Uq^ just equal to the wave function Uq in the absence of regulators. 
This contribution takes the form 

G'/'<Q{t) = {2nff_JhaHti)V''''{Tp{h))lm[ul{t)u;\h^ 

X / (fq |«q(ii)|' . (40) 

Jq<Q 

For q > Q the individual terms in Eq. (39) are quadratically divergent, 
so here we need an extended version of the WKB approximation (30), in 
which we keep terms in u of order and k^^I"^ as well as This is 

complicated by the presence of the potential term in Eq. (29), which couples 
wave functions with different ks. We will deal with this by considering the 
potential term in Eq. (29) as a perturbation. Of course, V"{}p) is not a per- 
turbation; it is of zeroth order in the loop-counting parameter g introduced 
at the end of Section II. However, each insertion of V"i^^ in the loop in 
Diagram II lowers its degree of divergence by two units, so the only terms 
we need to consider are those of zeroth and first order in V" which in 
the absence of cancelations are quadratically and logarithmically divergent, 
respectively. Terms of higher order in V" {Tp) are convergent even in the 
absence of cancelations, and are therefore negligible. 

To evaluate the terms in Gp^'^'^(t) of zeroth order in F"(^), we note 
that in the absence of the potential, u^{t\) is proportional to ?>kn'- 

^Nqiil) = ^NKUNq{t\) , (41) 

where 

UNq + SHilNq + {q^/a^)uNq + M^UNq = . (42) 

This contribution is 

G^/'>Qfi{t) = (27r)3£^diia3(ii)y""(^(ii))lm(^x2(i)^i;2(ii)) 

xJ2Z^^ [ d^qWNq{ti)f . (43) 

N ■' 
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The integrand is given by an asymptotic expansion derived in the Appendix. 
For both q^/a^iti) and much greater than both H'^ih) and H{ti), we 
have 



1 



2KiVga^(27r)" 



H + 2H^ {H + 3i?^)Mjv 

1 ~ O 



(44) 

where, as before, n%q{ti) = (^g/a(ti)^ The integral over q converges 

because Etv^t^^ = Eiv^jv^A^ = 0. The sum over N receives contri- 
butions from terms where xn is a regulator field Xn or the physical field 
Xo = V- Adding these contributions gives 

Gli>Q'\t) = 7r^*^dtia=^(ti)y""(^(ii))lm(u2(t)n;2(t,)) 

^ Z-'M^ In Mn + {H{h) + 2H\ti)) - ^ Z'^ lnM„) 



Q 



^-(i^(.,)+2^^(.,))ln(-| 



a2(t 

The regulator-dependent term arising from diagram II that are of first 
order in V"{lp) can be calculated by applying the rules of the "in-in" for- 
malism a diagram like that of diagram II, but with a V" insertion in the 
loop. This gives 

Gl^'>Q'\t) = -{27rf j'_^dha^{h)V""(lp{h)) f_^dt2 a^{t2)V" [lp{t2)) 



(45) 



-00 

X E 

KLMNM'N' 



7-1 7-1 



q>Q 



Sq Re ul{t)ul{h)ul(ti 



[e{h-t2)u^q{h)u^\{t2)u'},q{h)u'^\{t2) + l^2] I. (46) 



(This contribution is produced only by terms in which both interactions 
come from the time-ordered product in Eq. (13), or both from the anti- 
time-ordered product. As in the case of diagram I, the other terms make a 
negligible contribution to the part of the integral with q > Q.) The individ- 
ual terms in Eq. (46) are only logarithmically divergent, so we can evaluate 
this using the leading term (30) in the WKB approximation. Following the 
same limiting procedure as for diagram I, we find 

(5//,>Q,i(^) = vr r dh aHh) V""(lp{ti))v"(lp{h)) Im [ul{t)u;\h 
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i M^lnM„-M^lnM„ 

n m 



+ 2^Z-ilnM, + ln(-^ 



(47) 



Total IPI Amplitude 

The complete contribution of the two one-particle irreducible diagrams 
is given by the sum of the terms (35), (38), (40), (45), and (47): 

Gl^'it) = -2{27rf J'_^dh a^ih)V"'(lp{h)) j'^^dt2 a\t2)V"' {Tp{t2)) 



X Re <^ «p(t)u*(tl)u*(t2) / d q Uq{tl)u*g(t2)Uq>{tl)u*,{t2) 

+{27rf dh a^{h)V"'(Tp{ti)) dt2 a\t2)V"'{Tp{t2)) 
X |«p(t)pRe |tt*(ti)«p(i2) j d^qu*g{ti)u*g,{ti)uq(t2)ug/{t: 
+(27r)3|_*^a3(ti)y""((^(ti))lm{tx2(t)u;(ti)} d\\uq{h)\ 



+7r j dti a^{ti) 



V"'{lp{h)y + V""(lp{h))v"{lp{h)) Im {ul{t)u;\h)} 



+7r 



V- 7-1 / M^lnM„-M^lnM^ \ , .^^-i 

mn \ n -'"m / n 

£^dh a^{h)V""(lp{h))lm {ul{t)u;\ti)} 



\nMn + In 



Q 



a{ti 



+ {H{t,) + 2H\t,)) (^--Y. lnM„ - In ' 

To repeat, q' = |q — p|, and Q is any wave number for which Q'^/a^{t) 
is much larger than and V'^lp) and /a?{t) and much less than all 
regulator masses. In this range, the Q-dependence of the first and third 
terms is canceled by the explicit Q-dependence of the fourth and fifth terms. 

VI. CANCELING THE REGULATORS 

The terms in the counterterm (17) that are quadratic in the fluctuation 
make a contribution to the interaction-picture Hamiltonian of the form 



(48) 



^^quad^^) = lg(i) j d^x Sip\^,t) , 



(49) 
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where 



AV""{^) + 2B[V"'^{l^) + V" {lp)V"" {ip)] - 6CiH + 2H'')V"'{-^) 



(50) 

According to the rules of the "in-in" formahsm, this makes a contribution 
to the two-point function (33) given by 



AG'/'it) = 2(27r)3 r dh g{h)l^{ul{t)u;\h)} ■ 



(51) 



Comparing Eqs. (50) and (51) with (48), we see that in order to cancel 
the dependence of the one-particle irreducible two-point function on the 
regulator properties, we need 



A 
B 



1 



167r2 
1 



'Ml HMn/fiB) - Ml \li{Mm/iiB) 



(52) 



c = 



327r2 



M2 - M2 



(53) 
(54) 



(The first term in Eq. (52) does not depend on the units used for regulator 
masses in the logarithm, because Y^n^n^^n = 0-) Here /xa, fJ-B, and fic 
are unknown mass parameters. The presence of these parameters should 
not be seen as a drawback of this method; they reflect the real freedom we 
have to add finite regulator-independent terms to the original Lagrangian 
proportional to V"{(p) or V"'^{ip) or RV"{(p). 

Adding Eqs. (48) and (51) gives our final answer for the one-particle- 
irreducible part of the two-point function 

Gf^(t) + AGf^(t)= -2{2TTf f dha\h)V"'(ip{h)) r dt2a\t2)V" 

X Re < ul{t)u*p{h)u*{t2) / d^q Uq{h)u*q{t2)Uq' {h)u*q, {t2 
I Jq<Q 

+7r £ dh a^{h)V"'(ifih)yim{ul{t)u;\h)}ln (^^^^ 
+(27r)6 J^^ dh a'{h)V"'(ifih)) dt2 a^{t2)V"' (ip{t2)) 
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x\Up{t)\ Re<ul{ti)Upit2) / d'^qU*q{ti)u*q,{ti)Uq{t2)Ug'{t2) 



+ 



dh a'{h)V""(Tp{ti)) l^{ul{t)u;\h)}< 



'q<Q 



+ ^"(^^*^)) (^) - (^^'^^ + ^^^(^^)) (^) ^ ^^)^ 

For Q'^/a?{t) much larger than i?^, i/, and p'^/a?{t), all Q de- 

pendence cancels separately in the terms in square brackets on the first 
three lines and on the last three lines. In this form, the two-point function 
(including also the onc-particlc-rcduciblc contribution discussed in the fol- 
lowing section) can be calculated even if all we have for the wave functions 
Uq{t') is a numerical approximation. 

VII. ONE-PARTICLE-REDUCIBLE DIAGRAMS 

We now turn to the one-particle-reducible diagram III. In this diagram 
the two external lines come together in a three-field vertex, with the third 
line terminating cither in a three-field vertex to which is attached a scalar 
loop or a one-field vertex arising from the part of the one-loop counterterm 
(17) that is linear in 5ip. This part of the counterterm is 

Aifj^°(i) = Fit) J d^x 6(p{x,t) , (56) 

with F{t) given by 



7"= -a^ 



AV'"(lp) + 2BV"{lp)V'"{ip) - C{6H' + l2H)V"\lp) . (57) 



This diagram requires special treatment, because the line connecting the two 
vertices carries zero three-momentum. For this reason, here we will delay 
integrating over the difference x of the spatial coordinate of the two vertices. 
The full one-particle-reducible contribution to the two-point function (33) 
is then 



Gi^«(i) = 2(27r)3Re j_^dh a^{h)V"' (lp{ti))ul{h)u;\h) 

X r dt2l{t2) I d^x f-(T{5^(0,ti)(5(/?(x,t2)}>o + (Mx,ti)MO,t2)>o 
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where 

X{t2) ^ \a\t2)V"' irpit^)) f d'q Y: u§^{t2)u§r^{t2) + Ht2) . (59) 

KNN' 

In the first term in the square brackets in Eq. (58), both vertices come 
from the time-ordered product in Eq. (13), while in the second term, ver- 
tex 1 comes from the time-ordered product and vertex 2 from the anti- 
time-ordered product; in the complex conjugate time-ordered and anti-time- 
ordered products are interchanged. 

There is no problem here with ultraviolet divergences coming from the 
integral over q. Following the same procedure as in our treatment of diagram 
II in the preceeding two sections, we have 

+ 2JJ^te))ln(-^)+ri^ 

where Q is any wave number with Q'^/a?{t) much larger than H(t') and 
H^t') and V"(Tp{t')) for all t' < t. All dependence of Q cancels in this 
limit. 

But there is an apparent problem with infrared effects. Eq. (58) involves 
the integrals 

d'^x {T{6(p{0,h)S^{x,t2)})o and j d'^x {6<p{x,t2) S^{0,h))o ■ 

When we use Eq. (9) for the interaction-picture fields, the integrals over x 
pick out the value zero for the wave number q. But the wave function Ug{t) is 
not defined in the case q = 0, because in this case there is of course no time 
early enough so that q^ /a'^{t) is much larger than H'^{t) and \ V" {Tp{t)^\. For 
the same reason, the argument for the Bunch-Davies condition a(q)$o = 
breaks down for q = 0. 

Fortunately, we need the integrals over x only in the combination 

j d^x [-(r{M0,ti)Mx,i2)})o + (Mx,t2)M0,ii))o] 

= ie{ti-t2)G{tut2) (61) 

where 

G{ti,t2)^i j d'x ( [5(^(0, ti ), 599(x,t2)])^ . (62) 
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Despite the ambiguity in uo{t) and the inapphcability of the Bunch-Davies 
condition for q = 0, the function G{ti, ^2) is perfectly well-defined. It is the 
solution of the second-order differential equation 



V"(ip{h 



G{ti,t2) = 0, 



(63) 



subject to initial conditions dictated by the commutation relations (7) and 
(8): 

G{t2,t2) = 0, (64) 



tl=t2 



a-^{t2 



(65) 



The only property of the vacuum state used here is that it has zero momen- 
tum and unit norm. The general solution is 



G{ti,t2) = u{ti)u{t2) 



dt 



where u{t) is any solution of the q = wave equation 

u + ZHii + V"{lp)u = , 



(66) 



(67) 



that does not vanish between ti and t2- (For instance, for a general poten- 
tial and a de Sitter metric, we can take u = Tp, which does not vanish in 
typical inflationary models.) Putting this together, we have the one-particle- 
reducible contribution to the two-point function (33): 

rt 



G 



IPR 



-2(27r)^ dh a^h) V'"(lp{h)) lm{ul{t)ul\h)} 

dt2G{hM)At2) ■ (68) 

VIII. THE ONE-POINT FUNCTION 

In Section II we defined 5ip as the departure of the field (f from its 
classical value ^, not from its mean value, so we must expect (^95 to have a 
non-vanishing expectation value. As wc will see, this is closely related to 
quantities calculated in the previous section. 

According to the general diagrammatic rules, the vacuum expectation 
value of the Heisenberg picture scalar field fluctuation in one-loop order is 

(5(/'H(y,i))vAC°°^ ^ J ^^^1 / {S'Piy^t) S'P{^i,h))oI{h) + C.C. , 

(69) 
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with X given by Eq. (60) representing the insertion of a loop or a counterterm 
at the end of the single incoming line. In the term shown in Eq. (69) the 
single vertex comes from the time-ordered product in Eq. (13); in its complex 
conjugate, the vertex comes from the anti-time-ordered product. The two 
terms together involve the commutator of the field perturbations, so the 
one-point function may be written in terms of the function G defined by 
Eq. (62): 

{^^u{yA)TKc"^ = - r G(t,ii)x(ii) . (70) 

We see now that the contribution (68) of the one-particle-reducible diagrams 
to the two-point function may be simply expressed in terms of the mean 
fluctuation: 

Gf\t) = fjt,aHt,)V"'{Wi)) 

X (<5^H(0,ti))^Xc°°' lni{ulit)ufih)} . (71) 

This is the same as would be given by adding an interaction obtained by 
shifting dip by its expectation value: 

AHrit) = ^aHt)V"'(Tpit)) (5(^^(0, t))°Xc°°'' / d'x V(x,t) ■ (72) 
IX. INFRARED DIVERGENCES? 

Although the model treated in this paper is intended to provide an il- 
lustration of a method of dealing with ultraviolet divergences, it may be of 
some interest to look into the possible presence of infrared divergences in 
this model. For any fixed co-moving wave number q, the evolution of the 
wave function Uq{t) defined by Eqs. (11) and (12) becomes (^-independent 
once q/a{t) drops below H(t), so the behavior of the wave function for fixed 
t and g' — 7- is determined by the behavior of V"(Tp{t')^ and H{t') for t' — t- 0. 
We can distinguish two cases in which this problem is greatly simplified. 

Expansion-dominated : 

If V"(Tp{t')^ < if2(t') for t' 0, then as long as this inequality is satisfied, 
we can drop the potential term in Eq. (11), which then becomes the same 
as the differential equation for tensor fluctuations. It is well known[10] 

in this case that if H{t') -eH'^it') as t' 0, then the wave function 



Uq{ti) at a fixed time ti goes as q for q/a 

dependence is unaffected even if (t) drops below 



(ti) < H{ti). This q- 



V 



at some time 
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after q/a drops below H, since the evolution of the wave function at such 
times is g-indcpcndent. So (taking e < 1) the integral over q of the product 
Uq{ti)u*{t2) in the propagator will be infrared divergent if and only if e > 0. 
(We have been assuming that as time passes fluctuations leave the horizon 
rather than entering it, so this discussion is limited to the case e < 1. For 
the case e > 1, see ref. [11]-) There is no infrared divergence in the unlikely 
event that the expansion rate increases at very early times. 

Potential-dominated: 

If V"(Tp[t')j » for t' — > 0, then as long as this inequality is satisfied, 
Eqs. (11) and (12) have a WKB solution 

- (2.)3/^„3/4,)^ '''po/r-('">'«") . (73) 

where T is arbitrary, and 



Once q/a{t') falls below \ V" \ Tp{t')j |, the wave function Uq{t') becomes inde- 
pendent of g, aside from a g-dependent phase that is independent of t' . Later, 
H^{t') may or may not become comparable to or greater than V"(jp{t')j , 
but this cannot affect the g-dependence of the wave function. Therefore 
when the potential dominates at very early times, the product Ug{ti)u*{t2) 
in the propagator at fixed times ti and t2 becomes g-indcpcndent for q ^ 0, 
and there is no infrared divergence when we integrate the propagator over 

q- 

X. FURTHER ISSUES 

The method described here can of course be applied in this model to all 
one-loop correlation functions. The same counterterms, given by Eqs. (16) 
or (17) and (52)-(54) will remove dependence on the regulator properties, 
because the only ultraviolet divergences in one-loop one-particle-irreducible 
diagrams occur in the one-point and two-point functions, which we have 
already discussed in Sections V through VIII. The only ultraviolet diver- 
gences in higher correlation functions arise in diagrams in which trees are 
attached to loops at either one or two vertices, and the divergences in these 
loops are just those with which we have dealt. Multi-loop graphs are more 
challenging. 
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Beyond the simple model discussed here, of a scalar field in a fixed met- 
ric, there is the more realistic problem of scalar and tensor fluctuations in 
a theory of coupled scalar and gravitational fields. This is more compli- 
cated, because even in one-loop order there are quartic as well as quadratic 
and logarithmic ultraviolet divergences. That alone should not prevent the 
method described here from being applicable to realistic theories, at least 
for one-loop graphs, since divergences of any order can be eliminated by 
including enough regulator fields. 

A more serious problem is the difficulty of introducing regulator fields for 
the graviton propagator. (This problem is of course avoided in theories with 
large numbers of matter fields, where matter loops dominate over graviton 
loops.) If the only vertices that involve gravitons have a single graviton line 
attached to matter lines, then we can introduce regulators for the gravi- 
ton propagator by coupling heavy tensor fields with suitable Z-factors to 
the energy-momentum tensor. But it is not clear how to deal with graphs 
containing vertices to which are attached two or more graviton lines. 

This raises the question whether Pauli-Villars regularization is really 
necessary. The final results (55) and (68) for the one-particle irreducible 
and reducible parts of the two-point function could almost have been guessed 
without introducing regulator fields. It would only be necessary to intro- 
duce an ultraviolet cut-off at a sufficiently large co-moving wave number 
Q, calculate the Q-dependence of the resulting two-point function by using 
the WKB methods described in this paper, and then introduce a countert- 
erm of form (16), with A, B, and C chosen as functions of Q to cancel the 
Q-dependence found in this way. (This is not the adiabatic regularization 
procedure mentioned in Section IV, even though both procedures use WKB 
methods, because with a cut-off at Q only the part of the integrand for 
internal wave numbers larger than Q is affected.) Of course, this procedure 
leaves finite terms in A, B, and C undetermined, but they are undetermined 
anyway, since they represent the real possibility of changing the original La- 
grangian by adding corrections to the potential and adding a coupling of 
the scalar field to the spacetime curvature. The cut-off introduced in this 
way would not respect general covariance, but apparently one would get the 
correct results (55) and (68) anyway. 

There is something mysterious about this. The actual calculations in 
this paper were done for a fixed Robertson- Walker metric, Eq. (2). They 
would have been done in the same way by someone who had never heard 
of general covariance. Yet the infinities turned out to depend on H and 
H only in the combination H 2H^, proportional to the scalar spacetime 
curvature. We can understand this for a generally covariant regularization 
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procedure, like Pauli-Villars regularization, because in that case general 
covariance is broken only by the background, which presumably does not 
affect ultraviolet divergences. But how do these calculations know that they 
are supposed to give infinities that can be canceled by counterterms that 
are generally covariant, when we use a non-covariant cutoff on the internal 
wave number instead of introducing regulator fields? 
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APPENDIX: THE EXTENDED WKB APPROXIMATION 

Wc wish to find an asymptotic expression for the solution Uq{t) of the dif- 
ferential equation 

Uq{t) + 3H{t)Uq{t) + [q^/a'^{t))uq{t) + M^Uq{t) = (A.l) 

subject to the initial condition, that for t — > 0, 

(The effects of the potential arc treated separately in Section V.) We are 
interested in the behavior of Uq{t) at a fixed time t, when q/a{t) is much 
larger than H(t), but not necessarily greater than M. 
As an ansatz, we take 

(A.3) 

with /, g, etc. of zeroth order in q and M, and 

K{t) = ^q^/a?{t)+M'^ . (A.4) 

This clearly satisfies the initial condition (A.2). The differential equation 

(A.l) is satisfied by (A.3) to order k'^/^ and k^/^, while the terms in (A.l) 
of order n^^^"^ (counting M as being the same order as k) give 

d ff\ if- „ . 3i/2Af2 5M^H^ HM^\ 
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The terms in (A.l) of order k^^^"^ are more complicated, but fortunately 
we only need these terms in and for this purpose we can avoid having 

to work out these terms by using the time-dependence of the Wronskian: 

^ (A.6) 



Using (A. 3) gives 



^2id_ /Im/\ 

K dt \ K, J 



Uqllq 



.1/1 



-2i 



4iRe/ 



2V-^^ j^ + 0{k 



(A.7) 



Now, Eq. (A. 5) shows that d/dt(f/K) is imaginary, so since f{t)/K{t) 
vanishes for t — )• 0, f{t)/n{t) and hence f{t) is imaginary for all t. The 
first term on the right-hand side of Eq. (A.7) is constant, and the second 
term vanishes, so the constancy of this quantity requires the vanishing of 
the terms of order 



|/|2 + 2Re<7 



d flraf 
"'di 



But this is just what we need, for Eq. (A. 3) (with / imaginary) gives 



2K(t)(27r)3a3(t) 



1 + 



|/(t)p + 2Reff(t) 
^2(t) 



(A.9) 



Together with Eqs. (A. 5) and (A.8), this gives the result used in evaluating 
diagram II in Section V. 



1 



2«;a3(27r) 



2^2 
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